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WASHITA  RIVER  WATERSHED  RESEARCH  FINDINGS 


INTRODUCTION 

The  work  of  the  Southern  Plains  Watershed  Research  Center  is 
continuing  to  move  ahead  dynamically.  T3^ical  of  our  time,  1967  has 
been  a  year  of  change — a  year  that  has  produced  new  problems  and  fresh 
research  challenges,  like  the  water  quality  problems. 

This  report  highlights  our  research  activities  and  accomplishments. 
In  1961  the  Agricultural  Research  Service  started  a  hydrologic  research 
study  of  a  1,130-square-mile  area  in  a  central  segment  of  the  Washita 
River  Basin  in  Oklahoma  to  determine  the  change  in  dovmstream  runoff, 
sediment  flows,  and  groundwater  levels  when  flood  control  measures  are 
applied  to  the  uplands.  These  measures  include  land  treatment,  but  the 
dominant  feature  is  the  flood-detention  reservoir.  Detailed  rainfall, 
runoff,  sediment  yield,  and  groundwater  data  are  being  collected  and 
analyzed.  These  findings  are  of  great  interest  to  x^jater-short  areas 
such  as  Oklahoma. 

The  results  reported  here  are  based  on  detailed  field  observations 
being  made  by  the  Agricultural  Research  Service  in  cooperation  X'jith 
State  and  Federal  agencies  at  selected  locations  in  a  portion  of  the 
Washita  River  Basin  in  west-central  Oklahoma.  The  attached  figure 
shows  the  portion  of  the  Washita  River  Basin  selected  for  study. 


RAINFALL  PATTEPI^S  IN  THE  RESEARCH  AREA 
Arlin  D.  Nicks 

A  dense  rain  gage  network  has  been  Installed  on  the  1,130-sauare- 
mile  study  reach.  The  basic  network,  containing  168  recording  rain 
gages  spaced  on  a  3-  by  3-mile  square  grid  pattern,  provides  needed 
data  for  the  major  tributary  watersheds  and  the  entire  study  reach. 

The  grid  pattern  for  the  basic  network  is  shoxm  in  Figure  1.  Additional 
gages  have  been  installed  in  a  still  denser  pattern  on  smaller  water¬ 
sheds,  which  range  in  size  from  a  few  acres  to  7.6  square  miles. 

The  rain  gage  netx^ork  is  used  to  obtain  an  accurate  estimate  of 
area  rainfall  and  a  definition  of  storm  patterns  through  point 
measurements.  Occasionally,  the  runoff  response  from  a  storm  has 
indicated  that  a  rainfall  pattern  may  have  occurred  which  was  different 
from  that  defined  by  the  nccx^jork.  Additional  point  measurements 
between  stations  of  the  existing  netXrTork  could  have  defined  the 
correct  rainfall  patterns. 
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Figure  1 

The  rain  gage  distribution  over  the  study  reach 
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The  normal  annual  rainfall  for  the  region  according  to  Weather 
Bureau  records  (30-year  period,  1931-1960)  is  31  inches  per  year. 
However,  for  the  past  six  years  the  rainfall  has  been  below  normal. 
Table  1  shows  the  annual  rainfall  for  the  netX'Jork  for  the  period  1962 
to  1967: 


Tab le  1 ,  — Annual  P.ainf all 


Year 

Amount  (inches) 

1962 

28.44 

1963 

19.^1 

1964 

28.77 

1965 

26.40 

1966 

21.69 

1967 

26.41 

During  this  entire  72-month  period  only  18  months  have  had  rainfall 
above  normal. 

Weather  Bureau  records  indicate  that  annual  rainfall  varies  from 
from  35  to  23  inches  in  an  east  to  i/est  direction  across  the  study 
reach.  The  observations  during  the  past  six  years  show  that  there 
was  an  extreme  variation  of  6.32  inches  across  the  region;  however, 
the  minimum  rainfall  was  in  the  southwest  and  the  maxlmtum  in  the  north¬ 
east.  Figure  2  shows  the  average  annual  rainfall  over  the  region 
determined  from  the  measurements  obtained  from  the  rain  gage  network. 

The  major  portion  of  the  annual  precipitation  is  received  during 
the  period  of  Farch  to  October,  with  peaks  occurring  in  May  and 
September.  Both  convective  and  frontal  storms  have  produced  a  wide 
variation  in  rainfall  between  points  closely  spaced  on  the  network. 
Differences  as  large  as  4.5  inches  in  6  miles  have  been  observed  in 
both  local  and  widespread  storms^.  Yearly  occurrence  of  such  variations 
appears  to  be  characteristic  of  rainfall  for  this  climatic  region. 

Because  of  the  large  variation  in  rainfall  over  such  a  short 
distance,  a  study  was  made  to  determine  the  adequacy  of  a  network  for 
assessing  the  mean  daily  rainfall  over  the  entire  study  area.  The 
results  indicated  that  at  the  5-percent  level,  the  average  rainfall 
for  the  area  computed  from.  5  gages  was  significantly  different  from 
the  average  from  158  gages.  At  a  density  of  10  equally  spaced  gages, 
there  was  no  significant  difference,  indicating  that  daily  average 
rainfall  estimates  for  the  1,130-square-mile  area  could  be  obtained 
from  a  network  that  is  less  dense  than  that  installed. 


^Nicks,  A.  D.  and  Hartman,  Th  A.  Variations  of  rainfall  over  a  large 
gaged  area.  Transactions  ASAE,  Vol.  9,  No.  3,  pp  437,  439.  1966. 
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Figure  2 

Average  annual  rainfall  in  inches 
1962»1967 
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The  10“^age  network  might  supply  an  adequate  estimate  of  the 
daily  rainfall  for  the  1,130-square-mile  drainage  area,  but  for  small 
areas  in  the  basin  a  much  denser  network  is  required.  A  good  example 
of  this  need  is  the  storm  that  occurred  Vay  20,  1966  in  the  southwest 
corner  of  the  network  between  the  cities  of  Anadarko  and  Chickasha. 

The  network  observations  indicated  the  rainfall  pattern  shown  in  solid 
lines  on  Figure  3.  By  including  a  3.2-inch  rainfall  amount  reported 
by  a  resident  in  the  area,  the  rainfall  pattern  would  be  like  that 
shown  by  the  dotted  lines  in  Figure  3. 

Streamfloxi?  rises  sometimes  indicate  heavy  rainfalls  occurred  in 
the  watershed  that  were  not  recorded  on  any  of  the  gages  within  the 
network.  An  example  of  this  is  the  storm  of  September  3,  1966  in  the 
upper  West  Bitter  Creek  watershed.  A  5-foot  rise  occurred  at  the 
gaging  station  on  the  night  of  September  6,  yet  the  rainfall  data 
from  the  network  did  not  indicate  that  such  a  rise  was  possible.  So 
an  effort  was  made  to  collect  additional  rainfall  data  within  the  area. 
The  measurements  from  a  number  of  5-inch  glass  rain  gages  located  within 
the  network  were  added  to  those  of  the  network  gages.  The  map  in  Figure 
4  shows  the  resulting  rainfall  pattern  for  both  the  network  and  the 
glass  gage  data. 

Shox^m  in  Figure  5  is  the  storm  isohyetal  pattern  of  a  record 
rainfall  event  that  occurred  on  Sugar  Creek  watershed  in  the  northwest 
corner  of  the  network.  The  pattern  shown  v^as  produced  by  a  line  of 
thunderstorms  along  a  solid  front  that  moved  across  Oklahoma.  This 
storm,  which  follovred  two  days  of  near  uniform  rainfall  x^7ith  accumu¬ 
lations  up  to  4  inches,  caused  flash-flooding  and  considerable  property 
damage^ . 

A  noteworthy  characteristic  of  this  storm  was  the  high  rainfall 
intensity  recorded  at  various  stations.  Table  2  shows  the  maximum 
amount  recorded  for  time  intervals  of  5  minutes  to  6  hours.  In  most 
instances,  these  amounts  exceeded  the  maximum  recorded  for  the  state 
of  Oklahoma^.  The  clock-hour  distributions  are  shox'm  in  the  maps  on 
Figure  6. 


^Climatological  Data  National  Summary.  Weather  Bureau,  Vol.  16,  No.  9 
September  1965. 

3>taximum  Station  Precipitation  for  1,  2,  3,  6,  12,  and  24  hours.  Part 
]<XVI  Oklahoma,  Weather  Bureau  Technical  Paper  No.  15.  1961. 
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Figure  3 

Isohyetal  map  of  May  20,  1966  storm 
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Figure  4 

Isohyetal  map  of  September  3,  1966  storm 


Figure  5 

Storm  rainfall  map  for  Sugar  Creek  watershed 
September  20,  1965 
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Figure  6 

Hourly  rainfall  for  September  20j  1965  storm  on  Sugar  Creek  watershed 
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Table  2. — Maximum  intensities  for  selected  time  intervals  -  storm  of 


September  20,  1965 


Time 

Amount 

Intensity 

interval 

(inches) 

(inches/hour) 

5  min. 

.33 

9.96 

10  min. 

1.41 

8.46 

15  min. 

1.68 

6.72 

30  min. 

2,64 

5.28 

1  hr. 

4.61 

4.61 

2  hr. 

5,44 

2.72 

6  hr. 

8.83 

1.47 

Although 

the  rainfall  pattern  appears  to 

vary  uniformly  north-to- 

south  across 

the  watershed,  the  accumulations 

resulted  from  Intense 

cells  moving 

at  high  speeds.  Time-lapse  radar 

-scope  film  of  the  storm 

obtained  from  the  National  Severe  Storm  Laboratory  at  Norman,  Oklahoma^ 
showed  the  movement  of  small  cells  along  the  storm  line  to  the  north¬ 
east.  Figure  7,  a  photograph  of  the  PPI  scope  from  the  Norman  radar, 
shows  the  location  of  the  storm  with  respect  to  the  boundary  of  the 
network  at  2000  hours.  Tracking  the  storm  movement  by  tracing  con¬ 
secutive  frames  of  the  radar  film  reveals  an  average  cell  speed  of  30 
miles  per  hour  to  the  northeast.  A  similar  tracking  was  made  of  the 
surface  rainfall  by  plotting  average  10-minute  intensities  computed 
from  the  network  gages.  Maps  of  these  intensities  for  a  90-minute 
period  in  the  early  stage  of  the  storm  are  shoT^m  on  Figure  8.  Two 
cells  were  traced  by  these  methods,  the  first  occurring  from  1850  to 
1920  hours  near  the  top  of  the  watershed,  and  the  second  cell  traced 
was  from  1950  to  2020  hours  in  the  southern  part  of  the  watershed.  The 
tracks  shown  in  Figure  9  were  established  by  plotting  the  center  of 
these  cells  and  measuring  the  distance  and  direction  of  movement. 

Tracks  derived  from  the  radar  photographs  and  from  the  rain  gage 
network  agree  quite  well  in  both  speed  and  direction  of  movement.  In 
both  instances  the  average  speed  of  the  cells  was  30  miles  per  hour  to 
the  northeast. 


Figure  10  shews  a  plotting  of  rainfall  intensity  versus  time  for 
gages  along  the  line  of  the  cell  tracked  from  1950  to  2020  in  Figure  9. 
Gage  51  is  located  at  the  southwest  end  of  the  reach  and  gage  16  is  at 
the  northeast  end.  The  peaks  in  this  illustration  are  assumed  to 
represent  passage  of  cell  centers  over  the  gages.  Though  it  is  likely 
that  the  centers  of  the  cells  did  not  pass  directly  over  the  gages,  it 
is  believed  that  the  observed  intensities  are  satisfactory  for  tracing 
cell  movement. 


‘^Kessler,  Edwin.  Radar  measurement  for  the  assessment  of  areal  rainfall 
Review  and  Outlook.  Water  Resources  Research,  Vol.  2,  No.  3.  1966, 
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Figure  7 

PPI  scope  photograph  of  September  20,  1965  storm  at  2000  hours 
CFrom  National  Severe  Storms  WSR-S7  radar) 
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Figure  8 

lO-minute  intensity  (inches  per  hour)  plotting  of  storm  rainfall 

September  20,  1965 


s 

r.  ■■ 


■  _  .  '  'i'y  s. .,-  .* 

:v"  "■ 


•V  - 


» 


•  If 


'/  '*•*'  I  \  '  J  ‘  -  >  I  »  . 


/ 


Figure  9 

Cell  Cracks  from  radar  and 
rain  gage  network  data 
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Figure  10 

Rainfall  intensity  plots  for  selected 
gages  along  the  track  of  cells  of  the 
storm  of  September  20,  1965 
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Summary  and  Conclusions 


Some  small  isolated  storms  occurring  in  the  Southern  Plains 
climatic  zone  are  small  enough  to  pass  betXT^een  stations  of  a  3-  by 
3-mile  gage  spacing.  These  undetected  storms  can  cause  inconsisten¬ 
cies  in  the  rainfall-runoff  relationship  for  a  watershed.  The 
frequencies  of  occurrence  of  these  events  and  their  hydrologic 
significance  are  not  known  at  present.  Additional  data  are  needed 
so  that  better  analytical  techniques  and  procedures  for  projecting 
point  rainfall  to  areal  estimates  can  be  developed. 


WATERSHED  RUNOPi? 

Russell  R.  Schoof 

Instrumentation  and  Data  Acquisition 

Since  the  beginning  of  this  research  project  in  1961 5  mean  daily 
flows  for  about  100  station  years  of  record  from  20  main-stem  and 
tributary  gaging  stations  have  been  compiled.  The  present  gaging 
network  shown  in  Figure  11  includes  5  stations  on  the  Washita  River 
and  14  stations  on  tributaries. 

Performance  of  the  bubble  gage  xi^ater  level  recorders  has  been 
improved  by  installing  the  orifice  in  a  giant  sand  point^,  4  inches 
in  diameter  and  5  to  20  feet  long.  The  large  sand  point  provides 
ample  storage  for  sediment  that  enters  during  storm  flows  and  adequate 
anchorage  for  the  orifice.  A  self-cleaning  sand  point  can  be  construc¬ 
ted  from  an  open-bottomed,  short  length  of  pipe  fastened  to  a  bridge 
pier  at  an  elevation  such  that  the  stream  bed  scours  below  the  bottom 
of  the  pipe  during  a  rise. 

Discharge  is  measured  with  the  Price  current  meter  on  a  rod 
suspension  when  the  flox^  can  be  XA7aded,  and  on  a  cable  suspension  when 
the  floT.7  can  not  be  x-jaded .  Field  tests  indicate  that  approximately 
4,3  percent  less  discharge  is  measured  xijith  rod  suspension  than  with 
cable  suspension  when  the  same  meter  rating  table  is  used^.  Previous 
investigators  using  data  from  laboratory  calibrations  of  current  meters 
found  a  difference  of  2  percent  in  the  opposite  direction. 


^Hunt,  Charles  G. ;  Price,  Harold  E. ;  Schoof,  Russell  R.  ”,  and  Hartman, 
Monroe  A.  A  giant  sand  point  for  bubble  gage  orifice  housing.  USDA, 
ARS  41-129,  1966. 

^Schoof,  Russell  R.  Stream  gaging  precision  related  to  meter  suspension 
and  procedural  error.  Multilithed.  Master’s  thesis,  Oklahoma  State 
University,  May  1965. 


_.5_ 


/ 


ri 


•  c.U 


<  7  ’■*! 


j 


:jiT. 


M 


T  ~ 


•ii! . 

■:i  ‘  --  j .  v7  '■ 

:’  ■  •.  .  .‘t  •’  ~>c7rUL-‘  . 


i'i- 


r> 


{” 


*1  . 


i-'  ' 


^  •  ,•  ;  •  h 


.'•-hox 


-  i. 


>  y 


Figure  11 

Gaging  station  locations  on  the  Washita  River  study  reach 
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The  reasons  for  the  unexnected  result  are  not  known.  However, 
when  comparing  the  result  with  the  findings  o^  other  investigators 
note  that  in  this  experiment  the  tests  were  conducted  under  field 
conditions.  In  the  field  the  stream  gager  stands  in  the  water  x^hen 
making  a  rod  suspension  measurement,  whereas  in  the  laboratory  rating 
the  meter  is  suspended  from  a  tow  car  and  pulled  through  still  xixater. 
Even  though,  in  the  field,  the  stream  gager  stands  in  a  streamlined 
position  at  arm’s  length  from  the  wading  rod,  his  body  may  still 
retard  the  flow  in  the  vicinity  of  the  meter. 

A  system  for  computation,  checking,  and  routine  analysis  of 
stream,  flow  and  sediment  transport  data  has  been  developed^.  Economy 
is  maximized  and  human  error  is  minimized  by  using  electronic  data 
processing  equipment  including  a  chart  reader,  card  keypunch,  computers, 
and  card  sorter.  The  system  flow  chart,  shown  in  Figure  12,  includes 
12  Fortran  computer  programs  and  is  applicable  to  weir  control  and 
channel  control  gaging  stations  with  or  without  loop  ratings. 

Precipitation  and  runoff  have  been  below  normal  on  the  1,130- 
square-mile  study  area  since  the  initiation  of  the  research  nroject 
in  1961.  Punoff  from  the  study  area  has  been  less  than  1  inch  each 
year  since  1962.  The  normal,  long-term,  annual  runoff  for  the 
Chickasha  area  is  3  to  4  inches,  and  the  30-year  normal  precipitation 
is  31  inches  at  Chickasha.  Average  precipitation  on  the  study  area 
the  past  six  years  has  been  25,6  inches  per  year.  Annual  precipi¬ 
tation  and  runoff  for  the  period  of  record  of  each  watershed  is  shown 
in  Table  3. 

Flow  Losses  and  Abstractions 


To  supply  a  request  from  farmers  for  irrigation  water,  about  1,950 
acre-feet  of  water  were  released  by  the  Fort  Cobb  Conservancy  District 
from  Fort  Cobb  Reservoir,  about  44  channel-miles  upstream  from  Anadarko, 
during  August  1964.  This  release  into  the  Washita  River,  which  was 
nearly  dry,  provided  an  opportunity  to  obtain  data  on  channel  losses 
x-rithin  the  74  river-mile  study  reach  from  Anadarko  to  .Alex®. 

The  water  i-jas  released  at  150  cubic  feet  per  second  for  the  first 
two  days.  The  leading  edge  of  the  flow  moved  through  the  study  reach 
at  an  average  velocity  of  0.93  mile  per  hour.  The  average  loss  rates 


'Schoof,  Russell  R. ;  Seely,  Edward  R.  ;  Hicks,  A-rlin  D,  ;  Edens,  Carlton 
D. ^  and  Allen,  Paul  B.  Electronic  data  processing  system  for  stream 
flow.  Pending  publication. 

®Schoof,  Rvussell  R.  and  DeCoursey,  Donn  G,  Conveyance  of  irrigation 
vrater  in  a  natural  channel.  Published  in  the  Proceedings  of  the 
Second  Annual  American  Water  Resources  Conference,  November  20-22, 
1966. 
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within  subreaches  of  the  study  reach  varied  from  a  negative  loss 
(a  gain)  to  a  loss  as  high  as  4.39  acre-feet  per  mile  day.  The 
average  loss  rate  from  Fort  Cobb  to  Alex  vms  about  0.58  acre-foot 
per  mile  day.  Only  about  one-fourth  of  the  water  released  was 
diverted  from  the  channel  for  irrigation.  Host  of  the  remainder 
was  lost  by  seepage  or  evapotranspiration.  A  small  amount  was  in 
channel  storage  when  a  storm  broke  the  drought. 

Streamflow  losses  in  the  Anadarko  to  Verden  reach  of  the  Washita 
River  were  also  measured  during  natural  storm  flows.  Sixteen  runoff 
events  x^ith  peak  mean  daily  inflows  ranging  from  500  to  10,000  cfs 
showed  flow  losses  varying  from  5  to  77  percent  of  the  dally  inflow^. 
However,  some  of  the  water  lost  during  the  peak  returned  to  the  stream, 
during  the  recession.  Figure  13  portrays  three  of  the  runoff  events 
and  also  shows  the  accumulated  flow  loss  (abstraction)  or  gain  during 
a  six-month  period.  An  upward  trend  in  the  accumulated  abstraction  or 
gain  curve  indicates  flow  returning  to  the  stream  from  temporary 
storage  in  banks  and  flood  plain  soils, 

A  detailed  analysis  of  the  surface  water  budget  for  the  period 
September  20-30,  1965  for  the  study  reach  is  shovm  in  Table  4.  The 
average  loss  rate  of  27.9  acre-feet  per  mile  per  day  x^ras  much  higher 
than  the  loss  rate  observed  during  the  release  of  x-xater  from  Fort 
Cobb  because  the  flow  x-xas  much  deeper  in  the  stream  during  the 
natural  event  than  it  had  been  for  the  period  of  water  release.  High 
floxxs,  xxhlch  saturate  the  banks,  reduce  abstraction  from  succeeding 
flox^xs.  This  is  observed  in  three  runoff  events.  Also,  abstractions 
from  storm  runoff  events  xxere  greater  in  the  summer  months  x^xben  evap¬ 
oration  loss  xxas  great  and  the  antecedent  floxj  xxas  loxx.  Maximum 
gains  were  associated  with  the  winter  and  early  spring  flox-xs. 

Hydrologic  Influences  of  the  Sugar  Creek  Flood  Control  Program 

The  effects  of  an  upstream  flood  prevention  program,  on  the 
hydrology  of  Sugar  Creek,  a  240-SQuare-mile  tributary  to  the  Washita 
River,  were  studied^*^.  Figure  14  shows  hydrographs  for  three  storms 
occurring  after  the  installation  of  detention  reservoirs  and  the 
expansion  of  improved  land  use  measures  com.pared  with  hydrographs  for 
the  same  storms  estimated  for  the  pretreated  conditions.  The  storm  of 


^Schoof,  Russell  R. ;  Hartman,  Monroe  A.;  and  Hunt,  Charles  G.  Deter¬ 
mining  streamf lox<x  abstractions  from  antecedent  conditions.  Presented 
at  the  General  Assembly  of  the  International  Union  of  Geodesy  and 
Geophysics,  Berne,  Switzerland,  September  25-October  7,  1967,  and 
submitted  for  publication  in  lUGG  Transactions, 

^^Hartman,  Monroe  A.;  Ree,  William.  0.;  Schoof,  Russell  R. ;  and 
Blanchard,  Bruce  J.  Hydrologic  influences  of  a  flood  control 
program..  Published  in  the  Journal  of  the  Hydraulics  Division, 
Proceedings  of  the  ASCE,  Vol.  93,  No.  HY35  May  1967. 
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Figure  13 

Anadarko“Verden  reach  inflow  and  accumulated  abstractions 
and  gains  for  selected  events  in  1965  and  1966 
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Table  4. --Washita  surface  water  budget  (acre-feet)  for  the  period 


September  20-30,  1965 


Station 

Inflow  - 

Outflow  - 

Channa 1 
storage 

=  Loss 

River 
mi  les 

Loss 

acre-feet/mi/day 

Anadarko 

Sugar 

Verden 

05,096 

13,981 

83,323 

1,362 

53 

14,339 

21.8 

11.0 

39.7 

Verden 

Line 

Turnpike 

33,323 

54 

71,744 

1,936 

9,697 

26.9 

32.8 

Turnpike 

W.  Bitter 

E.  Bitter 
L.  Washita 
Winter 

Alex 

71,744 

123 

155 

246 

240 

66,892 

2,036 

3,530 

30.0 

10.7 

Anadarko- 

Alex 

99,395 

66,892 

5,437 

27,566 

89.7 

27.9 

■  \ I 


/ 


> 


c.  .'A' 


.  U- 


Figure  14 

Comparison  of  Sugar  Creek  hydrographs 


'M  i' 


September  1965  nearly  equaled  the  design  storm  and  provided  a  severe 
test  of  the  flood  control  v;orks.  Flood  peaks  Txere  reduced  at  least 
one-half  by  treatment.  The  recession  portions  of  the  af ter-treatm.ent 
hydrographs  were  lengthened  and  flattened.  The  character  of  the  re¬ 
cession  flow  was  influenced  by  the  number  of  detention  reservoirs 
discharging  into  the  stream  following  each  individual  storm.  The 
volume  of  water  infiltrated  was  reduced  on  the  flood  plain  and 
increased  in  the  detention  reservoirs. 


UNIT  SOURCE  watershed  RUNOFF  STUDIES 
Edd  D.  Rhoades 

Despite  the  limitations  inherent  in  relating  the  performance  of 
the  isolated  small  watershed  to  its  performance  as  a  unit  in  a  large 
complex  watershed,  the  small  experimental  watershed  is  still  a  useful 
tool  for  the  research  hydrologist.  He  selects  them  carefully  so  that 
they  are  representative  of  a  particular  condition.  By  manipulating 
the  cover  on  a  small  watershed,  he  can  study  the  effect  of  land  manage¬ 
ment  on  its  hydrologic  performance;  or  he  can  use  the  experimental 
watershed  to  evaluate  the  runoff-producing  characteristics  of  a  particu¬ 
lar  geology-soil-cover  complex. 

The  Southern  Plains  Watershed  Research  Center  has  selected  eight 
small  grassland  watersheds  and  eight  small  cropland  watersheds  for 
intensive  investigation  of  the  effect  of  land  use  on  hydrologic  per¬ 
formance.  Detailed  soil  surveys  and  topographic  surveys  i<reve  made  to 
define  their  physical  features.  Rainfall  is  measured  with  one  or  more 
recording  rain  gages  on  each  small  watershed;  runoff  is  measured  X'jlth 
V-notch  weirs  or  flumes,  or  by  observing  storage  changes  in  farm  ponds; 
and  soil  moisture  in  the  surface  4,0-foot  profile  is  measured  with  a 
neutron  soil  moisture  probe. 

Four  of  the  grassland  watersheds  (R-1,  R-2,  R-3,  and  R-4)  range 
in  size  from  18  to  25  acres  and  are  located  about  17  miles  northwest 
of  Chickasha,  Oklahoma,  on  a  permeable,  well-drained,  fine  sandy  loam. 
They  are  adjacent  to  each  other  and  are  all  \\Tithin  the  same  general 
pasture  area.  Portions  of  these  watersheds  were  cropped  several  years 
ago.  Although  there  has  been  considerable  erosion,  there  are  no  gullies 
or  active  erosion  on  these  watersheds  now.  Runoff  is  collected  in  a 
farm  pond  at  the  lower  side  of  each  watershed.  The  volume  of  runoff  is 
calculated  from  the  stage-volume  relationship  of  each  pond. 

The  other  four  grassland  watersheds  (R-5,  R-6,  R-7,  and  R-8) 
range  in  size  from  18  to  27  acres  and  are  located  about  15  miles  north¬ 
east  of  Chickasha,  Oklahom.a,  on  soils  that  are  loamy  to  silty  clay  loam 
with  good  surface  drainage,  but  poor  internal  drainage.  Runoff  is 
measured  through  V-notch  weirs.  Watersheds  R-5  and  R-6  are  native 
grassland  sites;  watersheds  R-7  and  R-8  were  formerly  cultivated,  but 
were  changed  to  grassland  after  they  had  eroded  badly. 
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The  eight  cropland  watersheds  are  located  in  the  Washita  River 
alluvium  and  range  in  size  from  13  to  44  acres.  Watersheds  C-l,  C-2, 
and  C-7  are  dryland  row  cropped,  primarily  in  cotton.  Watersheds  C-3 
and  C-4  are  in  irrigated  cotton;  C-5  and  C-6  are  in  wheat;  and  C-8 
in  two-year-old  alfalfa. 

Discussion 


Grassland  -  Table  5  shows  rainf all-runof f-sedim.ent  nroduction 
relations  for  the  rangeland  watershed  for  1967.  Runoff  from  water¬ 
sheds  R-1  through  R-4  has  been  low  during  the  5-year  period  of  record. 
The  5-year  (1953-1967)  average  annual  runoff  from  these  four  watersheds, 
R-1,  R-2,  R-l,  and  R-4,  was  .06,  .25,  .11,  and  .22  inches,  respectively. 
Precipitation  was  approximate!}^  equal  to  the  long-term  normal  of  31.60 
in  1967  and  below  normal  in  the  other  years  of  record.  Wo  high- 
intensity  storms  have  occurred  during  the  period  of  record. 

Runoff  amounts  from,  grassland  x-ratersheds  R-5  through  R-8  were 
much  greater  than  from  R-1  through  R-4  in  1967.  Watersheds  R-5  through 
R-8  V7ere  activated  during  1966;  therefore,  only  one  full  year  of 
records  is  available.  Precipitation  in  1967  on  these  x^atersheds  was 
belox7  the  long-term  normal,  x^ith  values  ranging  from  27.01  to  27.98 
inches.  Native  grassland  watersheds  R-5  and  R-6  reacted  similarly 
with  1,16  inches  and  1.23  inches  of  runoff,  respectively.  However, 
runoff  on  the  formerly  cultivated  watersheds  R-7  and  R-8  x^as  much 
larger,  3.89  inches  and  3.05  Inches,  respectively.  Severe  erosion 
during  and  subsequent  to  cropping  of  x^atersheds  R-7  and  R-8  appear  to 
be  one  reason  for  the  greater  runoff. 

The  severe  soil  loss  for  watersheds  R-7  and  R-8  reported  in  Table 
5  X7as  primarily  from  "inter-grass-clump"  areas,  indicating  that  these 
X7atersheds  had  not  yet  stabilized. 

Cropland  -  Table  5  also  shox-js  runoff  and  sediment  production  for 
1967  for  cropland  x^atersheds  C-l  through  C-8.  These  xratersheds  are 
bottomland,  V7ith  slopes  generally  less  than  0.1  percent.  Farming 
operations  on  all  watersheds  are  consistent  X7ith  normal  procedures 
x^ithin  the  area.  In  1967,  the  runoff  from  xratersheds  C-3  and  C-4 
in  irrigated  cotton  and  C-7  in  dryland  row  crops  was  2.50,  1.39,  and 
1.91  inches,  respectively.  These  watersheds  normally  have  a  large 
volume  of  runoff  xrhen  compared  with  the  other  xratersheds.  Some  of 
this  runoff  is  from  non- tilled,  irrigation  tailout  ditches  and  turn- 
roxrs.  Runoff  and  soil  loss  from  vratersheds  C-5  and  C-6,  soxm  to 
wheat  and  farmed  alike  were  1,32  and  1.19  inches  respectively. 

Runoff  and  soil  loss  from,  watersheds  C-l  (dryland  roxr  crop)  ,  C-2 
(dryland  roxr  crop)  ,  and  C-8  (dryland  alfalfa)  were  very  low  compared 
xrith  those  from  other  xratersheds , 

Rainfall-runoff-sediment  data  from  the  cropland  xratersheds  have 
been  collected  for  only  txro  years  and  no  definite  trends  have  been 
established. 
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Table  5. — Surrmary  of  rainfall,  runoff  and  sed?.iTient  loss  from  each  of  16 


unit 

source 

watersheds  *>  1967. 

Watershed 

Area 

Rainfall 

Runoff 

Sediment  loss 

No. 

(acres) 

(inches) 

(inches) 

(pounds) 

R-1 

17. 7S 

31.87 

.09 

i/ 

R-2 

24.03 

30.97 

.44 

1/ 

R-3 

25.84 

30.29 

.16 

1/ 

R-4 

18.12 

32.05 

.42 

1/ 

R-5 

23.72 

27.98 

1.16 

2,150 

R-6 

27.22 

27.94 

1.23 

3,300 

R-7 

19.19 

27.01 

3.89 

44,250 

R-8 

18.46 

27.44 

3.05 

57,350 

C-1 

17.83 

26.76 

0.08 

690 

C-2 

32.28 

24.85 

0.16 

1/ 

C-3 

44.26 

25.60 

2.50 

111,150 

C-4 

29.33 

25.55 

1.39 

36,100 

C-5 

12.75 

24.30 

1.32 

3,400 

C-6 

13.00 

24,66 

1. 19 

3,800 

C-7 

26.52 

24.38 

1.91 

31,800 

G-8 

27.23 

24.19 

0.57 

250 

_1/  Not  measured 
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Chickasha  Pumping  Sampler 


Sediment  yield  from  11  of  the  16  unit  source  T.jatersheds  was 
determined  during  1967  by  analysis  of  water  samples  hand  collected 
during  periods  of  runoff.  Since  runoff  often  starts  within  a  few 
minutes  after  rain  begins,  it  was  impossible  to  be  there  early 
enough  to  obtain  a  sufficient  number  of  samples  from  several  water¬ 
sheds  at  appropriate  times  to  accurately  calculate  sediment  quantities. 
To  obtain  a  continuous  curve  of  sediment  concentration,  frequent 
sampling  during  the  early  portion  of  a  runoff  event  is  required. 
Automatic  sampling  equipment  is  needed  to  do  this  economically.  The 
only  automatic  samplers  available  were  the  large  and  expensive  point 
samplers  in  use  at  larger  stations.  Technicians  at  Chickasha,  there¬ 
fore,  designed  and  built  a  smaller  and  less  expensive  automatic 
pumping  sampler  suitable  for  these  small  watersheds.  It  is  capable  of 
obtaining  a  sufficient  number  of  xi7ater  samples  to  adequately  show  the 
shape  of  a  sediment  concentration  curve  for  most  storm  events. 

This  new  sampler  is  powered  by  a  regular  12- volt  D.C.  automotive 
wet  cell  storage  battery.  A  mercury  switch  on  a  float  in  the  FW-1 
gage  well,  located  a  fex-j  feet  upstream  from  the  weir,  activates  the 
system  when  flox7  commences  over  the  x/eir.  A  12-hour  electric  clock 
pox-rers  a  notched  time  selector  disc  that  controls  the  frequency  of 
sampling.  The  present  sampler  is  constructed  to  obtain  samples  at 
10-minute  intervals  during  the  initial  period  of  runoff,  increasing 
by  small  intervals  of  time  to  60  m.lnutes  during  the  final  period  of 
sampling.  The  sampling  interval  may  be  changed  by  changing  the 
notched  time  interval  selector  disc,  A  rotating  circular  tray  holds 
28  pint  sample  bottles. 

The  electrical  system  consists  primarily  of  a  12-volt  battery, 
electric  clock,  a  four-unit  electric  timer,  and  a  system  of  micro 
switches  and  solenoids.  As  the  time  interval  selector  disc  rotates, 
the  arm  of  a  micro  switch  slides  into  one  of  the  notches  on  the  disc. 
This  turns  the  four-unit  electric  timer  on,  and  the  rotating  tray  is 
advanced  to  the  fill  position  under  the  diverter  discharge  spout.  A 
small  self-priming  electric  pump  then  starts  pumping  X'xater  from  the 
stream.  The  first  portion  of  the  xjater  that  is  pumped  purges  the 
lines  for  about  30  seconds,  then  a  solenoid  positions  the  diverter  to 
receive  the  sample.  After  the  sample  is  obtained,  the  diverter  is 
repositioned  and  pumping  stops,  A  short  tick  mark  is  m.ade  on  the 
Fl’J-1  water  level  chart  during  the  pumping  operation  by  use  of  a 
solenoid  that  moves  the  FW-1  float  tape.  The  sampling  operation  is 
complete.  A  second  sample  is  taken  when  the  unit  receives  another 
electrical  impulse  from  the  time  selector  disc.  The  operation  is  con¬ 
tinued  until  28  samples  are  taken,  or  the  mercury  sv/itch  breaks 
contact  with  the  clock  or  is  reset  by  a  technician,  A  safety  sx^itch 
on  the  time  selector  disc  prevents  the  refilling  of  bottles  already 
filled.  The  unit  functioned  satisfactorily  during  eight  runoff  events 
sampled  in  1967.  Similar  units  installed  in  all  unit  source  xixatersheds 
would  insure  obtaining  sufficient  data  to  calculate  a  fairly  accurate 
estimate  of  sediment  production. 
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SEDIMENT  TRANSPORT  MTO  CHANNEL  STUDIES 


Paul  B,  Allien  and  Norman  H.  Welch 

The  amounts  of  sediment  carried  by  the  river  and  its  tributaries 
are  measured  at  all  runoff  measuring  stations.  These  sediment  trans¬ 
port  data  and  the  channel  research  data  are  used  to  determine  changes 
caused  by  upstream  conservation  works  and  to  develop  equations  of 
channel  behavior  and  sediment  transport. 

Sediment  Yields 


1.  Short-term  yields;  Annual  yields  from  15  watersheds  are 
shown  in  Table  6.  Mean  daily  loads  and  incremental  transport  rates 
v/ere  used  to  determine  the  annual  yields.  The  data  in  Table  6  were 
computed  from  19,916  suspended  sediment  samples  taken  over  the  6-year 
period.  Of  these  samples,  2,130  v/ere  sieved  and  the  particle  size 
distribution  was  determined  with  a  visual-accumulation-tube  apparatus. 

Bed  material  gradation  analyses  were  made  on  705  bed  samples. 

2.  Long-term  yields;  No  analysis  has  been  made  of  long-term 
yields.  However,  the  long-term  yields  are  in  all  probability  greater 
than  the  1962-67  short-term  yields  shown  in  Table  6,  This  period 
includes  several  drought  years  and  lacks  large  out-of-bank  flox/s. 

Table  6. — Annual  sediment  transport  for  Washita  River  and  tributary  gaging 


stations ,  1962-67. 


Station 

Drainage 
area 
(sq .mi . ) 

Apprx . 
pet. 
sand 
in 

niGSs  • 

load 

Thousands 

of  tons 

1962 

1963 

1964  1 

1965 

1966 

1967 

Washita  River 

3,656 

9 

1,300 

172 

574 

951 

73.4 

92.5 

at  Anadarko 

Washita  River 

4,083 

18 

1,430 

172 

668 

1,120 

110 

230 

at  Verden 

Washita  River 

4,328 

13 

1,470* 

185* 

591 

1,200 

136 

284 

at  Chlckasha 

Washita  River 

4,707 

15 

1,050 

1,780 

265 

478 

at  Tabler 

Washita  River 

4,787 

17 

2,290 

348 

1,260 

2,030 

292 

534 

at  Alex 

Tonkavra  Creek 

26.0 

45 

1,7 

17.2 

4.6 

0.4 

2.3 

Sugar  Creek 

201.5 

67 

69.1 

257 

40,7 

124 

Delaware  Creek 

40,1 

39 

6,5 

4.2 

0.2 

3.2 

Line  Creek 

53.4 

3,3 

3.2 

2.5 

0.2 

11,9 

West  Bitter  Creek 

60.8 

9 

64.5 

13.6 

44.9 

43.3 

65.7 

East  Bitter  Creek 

35.0 

12 

65.1 

46,2 

39.8 

67,8 

Bedingfield  Creek 

19.2 

29.2 

25.2 

48.3 

Little  Washita  R. 

207.7 

30 

429 

172 

30.9 

65,5 

Big  Dry  Creek 

7.6 

1.2 

2.2 

0,2 

Winter  Creek 

33,3 

49 

141 

42,4 

7.5 

12,4 

*Data  from  old  station  5.7  miles  upstream. 
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3,  Effect  of  flood-detention  daros:  Seditnent  yields  were  only 
one-tenth  as  large  after  dams  were  installed  on  -Jinter  Creek.  This 
was  determined  xjith  a  double-mass  yield  curve  of  data  from  Winter  Creek 
and  the  adjacent,  nontreated  East  Bitter  Creek  watershed.  Such  a  large 
reduction  is  not  expected  in  the  long-term  yield  on  Winter  Creek 
because  the  analysis  period  was  generally  rainfall  deficient.  Reductions 
of  this  magnitude  are  not  expected  on  all  watersheds  because  the  Winter 
Creek  reduction  included  a  decrease  in  channel  erosion  as  well  as 
storage  of  upland  sediments  in  reservoirs. 

Mechanic s  of  S ediment  Transport 

1.  Sediment  and  velocity  distribution;  To  compare  total 
sediment  load  computations  made  with  the  Brooks,  modified  Einstein,  and 
Colby  procedures,  several  sediment  and  velocity  distribution  studies 
were  made  in  1963  and  1964  on  the  Washita  River.  Tables  7  and  8  give 
a  summary  of  the  sediment  and  flow  data  for  the  studies. 


Table  7. — Summary  of  sediment  and  velocity  distribution  data 


Washita 

River 

Station 

Disc. 

cfs 

Vel. 

fps. 

Ifeasured 
cone.-  ppm 

z 

von  Karman 
k 

.062-. 125  mm 

.125-. 250  mm 

Total 

Sand 

Range 

Av£. 

Range 

Avg . 

Range 

Avg . 

Anadarko 

1,630 

4.49 

5,580 

600 

.08-. 27 

.16 

.17-  .72 

.45 

.31-1.32 

.63 

Anadarko 

3,860 

4.53 

5,300 

820 

.08-. 26 

,14 

,25-  .69 

.40 

.55-1.33 

.95 

Verden 

4,360 

3.82 

4,890 

830 

.08-. 28 

.19 

.13-  .60 

.39 

.18-1.75 

.63 

Chickasha 

2,390 

3.48 

10,700 

590 

.03-. 20 

,14 

.17-  .44 

.35 

.36-1.07 

.60 

Alex 

1,250 

3.13 

8,800 

1,580 

.13-. 83 

.49 

.41-1.13 

.75 

.20-  ,88 

.43 

Alex 

4,250 

4,26 

7,490 

2,010 

.30-. 73 

.46 

.46-1.41 

.88 

.15-  .38 

,25 

Table  B.— Summary  of  total  load  data 


Sand  transport  rates 

Unmeas . 

cn 

tons /day 

load 

eo 

(by  Brooks) 

•H  G 

G 

QJ 

Pet. 

Pet. 

Washita 

C 

C 

Bed 

a) 

<4-1  4-1 

of 

of 

River 

Disc. 

Slope 

CO 

Fr. 

®50 

mmi 

■H  CO 

^  G 

sand 

total 

Station 

cfs 

ft/ft 

no . 

Me  as . 

O  ‘H 

S  W 

Colby 

Brooks 

load 

load 

Anadarko 

1,630 

0.0022* 

.039 

.40 

** 

2,460 

— 

3,200 

2,770 

11.2 

1,3 

Anadarko 

3,860 

0.00095 

.041 

.28 

** 

7,940 

— 

4,250 

8,100 

2.0 

0.3 

Verden 

4,360 

0.00036 

.032 

.22 

0.200 

10,400 

12,200 

13,800 

10,500 

1.0 

0.2 

Chickasha 

2,390 

0.00036 

.033 

.22 

0.180 

3,860 

4,620 

5,570 

3,970 

2.8 

0.3 

Alex 

1,250 

0.00032 

,025 

.24 

0.140 

5,290 

6,310 

7,750 

7,140 

26,0 

5.9 

Alex 

4,250 

0.00032 

.021 

.29 

0.140 

23,100 

27,700 

30,500 

29.300 

21.1 

^Approximate  slope  for  this  rock  control  reach 
**Solid  sandstone  bed 
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2.  Comparison  of  total  load  procedures:  Total  sand  transport 
rates  computed  by  the  modified  Einstein,  Colby,  and  Brooks  procedures 
are  shown  in  Table  8.  The  spatial  distribution  of  silt  and  clay  in 
the  measured  zone  of  the  stream  was  quite  good.  For  example,  the  silt 
and  clay  concentration  varied  only  from  about  4,200  ppm  near  the  sur¬ 
face  to  4,400  ppm  0.3  foot  from  the  bed  for  the  3,860-cfs  flow  at  the 
Anadarko  station.  From  this  it  was  assumed  that  the  silt  and  clay 
concentration  in  the  unmeasured  zone  was  essentially  the  same  as  in 
the  measured  zone.  Therefore,  the  measured  load  was  used  as  the  total 
load.  Because  the  Brooks  procedure  uses  field-determined  z  and  k  values, 
it  is  believed  to  provide  the  m.ost  reliable  estimate  even  though  it  does 
not  account  for  bed  load.  Calculations  indicate  that  bed  load  consti¬ 
tutes  less  than  2  percent  of  the  total  sand  load  or  less  than  1/2 
percent  of  the  total  load  for  the  studied  flows.  This  small  percentage 
is  due  to  the  fine  bed  material,  low  velocities  near  the  bed,  and 
relatively  narrow  channel  bed  widths.  For  example,  for  the  4,250-cfs 
flow  at  the  Alex  station,  if  the  thickness  of  the  moving  bed  is  taken 

as  twice  the  mean  particle  size,  the  velocity  of  the  bed  movement 
assumed  to  be  1  fps,  the  width  to  be  100  feet,  and  the  bed  material 
volume  weight  to  be  100  pounds  per  cubic  foot,  the  bed  load  would  be 
0.7  percent  of  the  total  sand  transport.  Table  8  shovrs  that  the 
Washita  River  in  the  vicinity  of  Verden  and  Chickasha,  Oklahoma  carries 
virtually  all  transported  material  as  suspended  load.  I-^^ere  the 
channel  is  wider,  as  near  Alex,  however,  roughly  25  percent  of  the 
sand  load  or  6  percent  of  the  total  load  cannot  be  measured.  It 
appears  from  the  Table  8  data  that  both  the  Colby  mean  velocity  method 
and  modified  Einstein  total  load  procedures  are  satisfactory  for  com¬ 
puting  total  load  in  the  lower  reach.  Because  each  Brooks  computation 
requires  lengthy  field,  laboratory,  and  office  work,  the  Brooks  method 
is  used  only  occasionally  to  check  the  accuracy  of  other  methods. 

3.  Effect  of  flow  temperature  on  transport:  For  3  days  during 
the  September  1965  bank-full  flow  on  the  Washita  Fiver  near  Alex,  a 
positive  correlation  existed  between  flow  tem.perature  and  the  fine 
material  in  transport.  This  was  believed  to  be  the  result  of  in¬ 
creased  bank  erosion  caused  by  an  increase  in  turbulence  with  warm 
water.  This  positive  correlation  between  temperature  and  transport 

of  fine  particles  is  especially  interesting,  because  other  researchers 
have  found  numerous  instances  of  negative  correlation  between  tem¬ 
perature  and  the  transport  of  coarse  material. 

Computational  Procedures 


1.  Continuous  sediment  concentration  method;  The  data  in  Table 
6  were  computed  by  taking  sufficient  suspended  sediment  samples 
during  runoff  events  to  construct  continuous  concentration  curves  on 
the  stage  recorder  charts.  Measured  transport  was  determined  by  a 
chart  reader  and  computer  using  the  relation 
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where 

T  =  transport  in  tons 

k  =  units  conversion  constant  =  0.0001123 

C  =  concentration  in  ppm 

Q  =  discharge  in  cfs 

t  =  time  in  hours 

2.  Unmeasured  load;  For  each  gaging  station  with  significant 
unmeasured  load,  the  modified  Einstein  total  load  procedure  was  used 
to  develop  a  discharge  vs.  unmeasured  load  curve.  Unmeasured  load 
determined  from  these  curves  was  added  to  the  measured  load  to  obtain 
the  total  load  data  shown  in  Table  6. 

3.  Sediment  rating  curve  vs.  continuous  concentration  curve 
methods;  A  comparison  is  being  made  of  sediment  discharge  from  small 
watersheds  (30  to  50  sq.  mi.)  as  computed  by  these  two  methods. 
Agreement  of  the  estimates  by  the  two  methods  was  erratic  when  only 
one  year  of  flow  and  sediment  data  was  used.  Differences  ranged  from 
8  to  100  percent.  Agreement  was  very  good,  however,  when  two  years  of 
flot7  and  sediment  data  were  used.  Differences  then  ranged  from  1  to 

4  percent.  Good  sample  coverage  is  especially  important  during  the 
higher  stapes  when  most  of  the  transport  takes  place.  Essentially  the 
same  flow  duration  curve  resulted  when  mean  daily  flows  were  used  with 
the  instantaneous  flow  peak  plotted  at  the  0.01  percent  tim.e  point  as 
when  time- incremental  (do^m  to  6  minutes)  flow  data  were  used. 

Accurate  definition  of  flow  duration  curves  is  needed  at  the  high-flow 
end  because  over  half  the  sediment  can  be  transported  between  the 
one-day  frequency  (0.27  percent  frequency)  and  the  0.01  percent 
frequency . 

Channel  Bed  Studies 

1.  Particle-size  distribution  of  bed  vs,  suspended  material: 
Because  the  gradation  of  the  bed  material  is  one  parameter  used  in 
computing  unmeasured  load,  an  understanding  of  bed  material  changes 
is  desirable.  In  1964,  bed  samples  were  taken  throughout  a  runoff 
event  at  Verden  that  peaked  at  6,300  cfs.  The  gradation  of  samples 
taken  with  a  BM-54  sampler  tended  to  follow  the  gradation  of  the 
suspended  material,  suggesting  a  continual  exchange  of  particles 
between  the  bed  and  the  suspended  load.  The  bed  contained  more  silt 
and  clay  v7hen  the  flow  was  near  the  peak  because  these  particles 
would  get  trapped  in  the  churning  bed.  During  and  after  the  recession, 
these  silt  and  clay  particles  V7ere  gradually  washed  from  the  bed. 
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2.  Gradation  of  bed  with  depth;  Short  reaches  of  the  Washita 
River  bed  have  been  sampled  intensively  in  cross  section  and  depth 
(5 ’“8')  at  Verden,  Chickasha,  and  Pauls  Valley.  The  median  oarticle 
size  varied  in  each  section  from  less  than  0.2  mm  near  the  surface 
to  greater  than  0.5  mjn  4  to  8  feet  below  the  surface  in  one  or  m.ore 
areas.  The  cause  of  these  coarse  areas  is  not  understood,  but  sub-bed 
flow  is  suspected  to  influence  their  development, 

Channel  Behavior  Studies 

1.  Aggradation  or  degradation;  During  the  past  6  years  numerous 
flow  events  with  associated  channel  scour  or  fill,  or  flood  plain  scour 
or  fill  have  been  observed  or  recorded.  Generally,  channels  and  flood 
plains  have  changed  little  in  watersheds  V7ith  principally  shale  geology. 
However,  large  channel  and  flood  plain  changes  have  occurred  in  sand¬ 
stone  geology  watersheds.  Intense  spring  thunderstorms  that  nroduce 
very  high  sediment  concentrations  and  less  than  bank-full  flows  have 
caused  the  most  channel  aggradation.  Slightly  out-of-bank  flows  on 
these  sandy  tributaries  cause  considerable  over-bank  deposition.  The 
extremely  large  8,600-cfs  flow  on  Sugar  Creek  in  September  1965  caused 
channel  bed  scour,  large  sand  splays  along  the  banks,  and  flood  plain- 
scour  at  various  places  in  the  valley. 

2.  Channel  shape  vs.  material  in  transport;  Streams  with  a  high 
sand  load  are  wide  and  shallow.  The  width- to-depth  ratio  (w/d)  is  de¬ 
fined  by  the  equation,  w/d  =  0.17X  +  3.55  where  X  is  the  percent  sand 
in  the  measured  load. 

3.  Sediment  progression;  The  intensive  sample  collection  during 
storm  flox7  at  the  sampling  stations  has  expanded  our  understanding  of 
the  progression  of  sediment  through  the  basin.  Continuous  sediment 
concentration  curves  of  large  flows  that  originate  entirely  above 
Anadarko  retain  almost  the  identical  shape  through  the  stable  channel 
past  Chickasha.  Below  Chickasha,  where  the  Washita  River  channel  is 
wider  and  less  stable,  there  is  a  fairly  large  pickup  in  sediment 
during  large  floxirs.  This  is  caused  by  the  flushing  action  of  the 
large  flows  on  sediment  stored  in  the  channel.  The  stored  sediment 
results  x^hen  large  tributary  flows  enter  the  Washita  River  when  the 
river  is  at  low  flow  and  deposit  much  of  their  load.  Cross  sections 
surveyed  by  the  Corps  of  Engineers,  U.  S.  Bureau  of  Reclamation,  and 
Agricultural  Research  Service  have  shovm  that  in  the  Tabler-Alex 
reach  of  the  Washita  River  the  bed  level  has  fluctuated  as  much  as 

o  feet  in  a  short  period  of  years  because  of  the  storing  and  flushing 
actions . 

4.  Bank  erosion;  Evidence  is  accumulating  that  excessive  sand 
loads  contribute  to  accelerated  bank  erosion  on  the  Washita  River. 
Although  the  mechanics  of  this  erosion  are  not  understood,  consider¬ 
ation  might  still  be  given  to  sand  transport  regulation  in  a  major 
bank-protection  plan. 
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5.  Channel  shape  vs.  vegetation;  A  vegetation  species  and 
population  count  was  made  at  several  points  on  the  Washita  River  to 
see  if  the  type  of  vegetation  affects  channel  shapes.  The  results 
were  inconclusive;  however,  a  plant  succession  was  observed.  As 
alluvial  deposits  age  and  build  higher  above  the  water  table,  the 
woody  plant  sequence  is  salt  cedar,  willow-cottonwood,  cottonwood, 
and  cottonwood-elra-hackberry . 

6.  Cross  sections;  To  document  channel  changes,  approximately 
200  channel  cross  sections  are  being  surveyed  periodically  on  the 
Washita  River  and  tributaries. 


Rimping  Samplers 

1.  XPS-62  pumping- type  sediment  sampler;  One  sampler  of  this 
type  has  been  used  for  two  years  at  two  Washita  River  stations.  After 
initial  m.alf unctions  x<rere  corrected,  the  sampler  has  operated  correctly 
and  has  proved  to  be  mechanically,  hydraulically,  and  electrically 
dependable.  The  greatest  problem  remaining  is  finding  a  location  for 
the  intake  so  that  pumped  sample^concentraticn  will  have  a  consistent 
relation  with  the  average  suspended  concentration  of  the  stream.. 

2.  PS-66  numping-type  sediment  sampler;  Two  of  these  were 
installed  on  the  Washita  River  in  1967.  Neither  collected  representa¬ 
tive  samples  with  the  original  2-minute  pumping  cycle  timer.  No 
samples  have  been  collected  since  the  2-minute  timers  were  replaced 
vjith  5-minute  timers. 

3.  Chickasha  pumping  sampler;  These  28-bottle,  small,  fully 
automatic  samplers  are  built  locally  for  small  watersheds.  They  are 
described  under  unit  source  watershed  runoff  studies. 


GEOMORPHOLOGIC  STUDY  OF  THE  WASHITA  RIVER  BASIN 
Don  W.  Goss,  Alex  P..  Ross,  and  Tommy  3.  Thompson 

A  knowledge  of  the  past  history  of  the  Washita  River  is  very 
important  in  understanding  and  predicting  future  trends.  To  obtain 
the  required  geomorphologic  history  and  geomorphic  form  of  the  river, 
several  lines  of  study  have  been  undertaken  recently.  Preliminary 
results  of  these  studies  will  be  discussed  briefly. 

To  determine  changes  in  alignment  and  length  of  the  Washita  River 
channel,  a  map  was  compiled  from  aerial  photographs  showing  the 
location  of  the  channel  in  1937,  1948,  1955,  and  1961.  In  addition, 
the  channel  locations  obtained  from  original  land  surveys  viere  added 
to  the  map.  Channel  lengths  for  each  year  for  which  data  are  available 
have  been  measured  on  the  map  and  the  results  suimnarizcd  in  Table  9. 

To  date,  no  maps  have  been  located  showing  the  Washita  Rivet  channel 
accurately  in  the  interval  between  the  original  1873  survey  and  the 
earliest  aerial  photographs  of  1937. 
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Table  9, — Channel  length  changes,  Washita  River,  in73~1961 


Year 

A 

B 

C 

1873 

42.3 

1898 

— — 

23.7 

14.0 

1937 

37.9 

25.2 

18.5 

1948 

38.0 

25.4 

17.5 

1955 

35.3 

25.6 

14.1 

1961 

36.5 

25.7 

13.6 

Column  A  Is  distance  from  Anadarko  gage  site  to  a  point  on  the  river  2, 
miles  west  of  Chickasha.  This  is  the  area  covered  in  the  1873  land 
survey. 


Column  B  is  distance  from  end  point  of  1873  survey  to  the  mouth  of  the 
Little  Washita  River. 

Column  C  is  distance  from  the  mouth  of  the  Little  Washita  River  to  the 
Alex  gage  site. 

Dating  of  the  alluvial  deposits  is  very  important  in  the  develop¬ 
ment  of  the  geomorphic  history  of  the  valley.  Animal  bones  recovered 
from  flood-plain  and  terrace  deposits  have  been  identified  by  personnel 
of  the  Museum  of  the  Great  Plains  at  Lawton,  Oklahoma  as  modern  cattle, 
canine,  and  bison.  It  was  evident  from  the  manner  of  occurrence  of  the 
bones  that  they  represented  burial  in-situ  shortly  after  death.  A 
group  of  bison  bones  sent  to  the  USDA  Sedimentation  Laboratory,  Oxford, 
Mississippi  for  carbon  dating  were  less  than  200  years  old. 

Two  buried  soils  from  lower  terrace  deposits  of  the  Washita  River 
V7ere  dated  by  L.  L.  TicDoxixell  at  the  USDA  Laboratory  at  Oxford,  Missis¬ 
sippi.  The  deepest  soil,  10  feet  belcni?  the  terrace  surface,  was  dated 
as  1,760  ±  150  years  Before  Present  (B.P.).  The  shallowest  soil,  8 
feet  below  the  terrace  surface,  was  dated  as  1,000  ±  100  years  B.P. 

Charcoal  samples  collected  7  and  9  feet  below  the  surface  of  the 
alluvium  were  dated  as  2,025  ±  150  and  2,215  ±  150  years  B.P.  respec¬ 
tively.  These  samples  V7ere  collected  in  the  Winter  Creek  channel  bank, 
about  3  miles  from  its  confluence  with  the  Washita  ’^iver. 

A  sample  of  wood  was  recently  recovered  from  a  depth  of  about  50 
feet  during  drilling  operations  near  Anadarko.  Samples  of  this  v70od 
were  dated  by  Isotopes,  Incorporated  as  3,850  ±  110  years  B.P. 

Mr.  A.  D.  Bull,  Soil  Conscx vation  Service  (retired),  CbleVasha, 
collected  some  walnut  from  Sandstone  Creek  (B.oge.r  Mills  County)  at  a 
depth  of  13.7  feet  and  some  Bois  d'Arc  from  Ctiner  (I'^oClaln  County) 

at  a  depth  of  8  feet.  These  samp lee  were  dated  as  1,025  ±  100  and 
1,548  ±  75  years  B.P.  respectively,  Lot-h  Samlstoue  r.reek  and  Criuex 
Creek  are  tributaries  of  the  Washita  River . 
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A  maminoth  has  been  uncovered  in  the  headwaters  of  Tonkavra  Creek o 
The  excavation  of  the  maminoth  was  directed  by  the  Huseum  of  the  Great 
Plains,  Lawton,  Oklahoma,  and  has  been  termed  the  Domebo  site.  From 
radiocarbon  datings  of  various  organic  remains  at  the  Domebo  site,  the 
age  of  the  site  was  concluded  to  be  about  11,200  years  B.P.  Table  10 
shows  the  known  radiocarbon  dates  of  materials  from  the  Washita  River 
Basin. 


Table  10. — Radiocarbon  dates  from  the  Washita  River  Basin 


Material 

Location 

Depth 

in  feet 

Age 

Years  B. 

P. 

Bison  bones 

Washita  River  bank 

14 

Recent 

Bison  bones 

Washita  River  bank 

4 

Recent 

Bison  bones 

Washita  River  bank 

— 

<200 

Buried  soil 

NT'J  1/4  SE  1/4  Sec.  .20, 

8 

1,000  ± 

100 

Buried  soil 

T.  6N.,R,  6W. ,  Grady  Co. 

MW  1/4  SE  1/4  Sec.  20, 

10 

1,760  ± 

150 

Charcoal 

T.  6N, ,R.  ,  Grady  Co, 

SE  1/4  ME  1/4  Sec.  25 

7 

2,025  ± 

150 

Charcoal 

T.  6N. ,R,  6W, ,  Grady  Co. 

SE  1/4  ME  1/4  Sec.  25 

9 

2,215  ± 

150 

Walnut 

T.  6F.,R.  6W. ,  Grady  Co. 
Sandstone  Creek 

13,7 

1,025  ± 

100 

Bois  d'Arc 

Criner  Creek 

8 

1,548  ± 

75 

Mammoth 

Tonkax'ja  Creek 

Streamhed 

11,200 

(Domebo  site) 
Wood 

Terrace  near  Anadarko 

50 

3,850  ± 

110 

Two  projectile  points  have  been  uncovered  from  alluvial  deposits 
in  the  Washita  River  terrace.  The  first  was  discovered  in  the  bank  of 
the  Washita  River,  and  Identified  by  Dr.  R.  E.  Bell  of  the  University 
of  Oklahoma  Museum  as  a  Washita  type  in  use  from  1400-1500  A.D.  Un¬ 
fortunately,  it  cannot  be  shown  that  this  represents  the  original 
burial  of  the  point.  The  second  projectile  point  was  uncovered  during 
coring  operations  at  a  depth  of  73  feet.  The  age  of  this  point  was 
placed  by  Dr.  Bell  at  about  7,000  years. 

In  order  to  describe  the  physical  characteristics  of  the  Spring 
Creek  v/atershed,  a  number  of  geomorphic  parameters  have  been  determined. 
This  x^ras  in  conjunction  with  a  subdivided  watershed  study.  After  run¬ 
off  data  are  obtained  from  the  xiiatersheds ,  the  relationship,  if  any, 
of  these  geopiorphlc  parameters  to  hydrologic  performance  X7ill  be 
determined.  The  parameters  chosen  are  those  that  have  been  accepted  by 
geomorphologists  to  be  definitive  or  indicative  of  some  character  of 
the  watershed.  These  parameters  are  area,  total  stream  miles,  drainage 
density,  stream  orders,  number  and  length  of  streams  in  a  given  order, 
basin  perimeter,  and  maximum  basin  length  parallel  to  the  major  drain¬ 
age  line.  Also,  the  circularity  ratio  and  the  elongation  ratio,  x-jhich 
are  dimensionless  parameters,  have  been  calculated.  The  parameters  for 
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spring  Creek  watershed  are  presented  in  Table  11.  Most  of  these 
parameters  have  also  been  determined  for  the  basins  of  stream  order 
three  or  greater  in  the  Little  Washita  River  watershed  above  Charlie 
Creek,  These  are  presented  in  Table  12. 

Dimensionless  parameters  are  quite  useful  in  defining  watersheds. 
Many  times  they  become  valuable  tools  in  the  development  of  mathematical 
relations  pertaining  to  hydrologic  characteristics  of  the  basin.  New 
dimensionless  parameters  have  been  Investigated  and,  thus  far,  two  have 
been  developed  that  warrant  additional  study.  They  have  been  defined 
as  the  length-width  ratio  and  the  basin  eccentricity.  Both  of  these 
parameters  are  derived  from  an  equivalent  elliptical  form  of  the  basin. 
Only  the  area  (A)  and  the  maximum  length  of  the  basin  parallel  to  the 
main  drainage  line  (L)  are  required  to  determine  this  equivalent  form. 

An  ellipse  is  defined  by  the  equation 

y2  y2 

2L  +  JL=:  1 

aZ  b2 


In  this  relation  ''a"  is  one-half  the  length  of  the  major  axis  and  ’’b" 
is  one-half  the  length  of  the  minor  axis  of  the  ellipse.  The  perimeter 
of  the  ellipse  is  defined  by  "X''  and  ”Y"  values  satisfying  the  equation. 
The  dimensions  of  the  equivalent  ellipse  for  the  basin  are  determined 


by  assigning  to 


the  value  —  and  by  calculating  "b**  from  the  relation 


A  =  nab,  where  A  is  the  known  area  of  the  basin.  The  ellipses  con¬ 
structed  to  scale  of  four  selected  third-order  basins  of  the  Little 
Washita  drainage  system  are  represented  in  Figure  15.  Figures  16  and 
17  show  the  comparison  of  the  elliptical  form  to  the  actual  shape  of 
each  basin.  It  can  be  seen  that  the  elliptical  form  is  very  repre¬ 
sentative  of  the  true  basin  shape. 


The  shape  of  an  ellipse  can  be  defined  in  a  comparative  form  by 
using  the  two  previously  mentioned  dimensionless  parameters.  The 
length-width  ratio  (W)  is  the  length  of  the  major  axis  divided  by  the 
length  of  the  minor  axis. 


This  value  will  be  one  for  a  circle,  and  approach  infinity  as  the 
ellipse  becomes  so  narrow  that  it  appears  to  be  a  straight  line.  Hio 
shape  of  an  ellipse  is  also  dependent  on  its  eccentricity  (E) . 


1 
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Table  11 . “-Spring  Crsek  watershed  geomorphlc  parameters 
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Table  12. — Geomorphic  parameters  for  certain  basins  in  the  Little  Washita 
River  watershed  above  Charlie  Creek 
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Table  12. --Gaomorphic  parameters  for  certain  basins  in  the  Little  Washita 


Rivar  watarshed  above  Charlia  Creak  (cont’d) 


Basin 

identification 

number 

Area 

(miles2) 

Number  of 
in  each 

1st  :2nd; 

L ? : 

streams 
order  1/ 

3rdj4th 

Circularity 

ratio 

Elongation 

ratio 

Basin 

perimeter 

(miles) 

maximum 

basin 

length 

(miles) 

Mean 

bifurcation  j 

ratio  1 

2.262 

10 

3 

1 

1 

.33281 .7529 

1 

I 

6.402! 

j 

2.254 

3.255 

^is'^e^i 

1.487 

11 

3 

1 

.  7409 

.5526 

5.8341 

] 

2.490 

3.513 

S5l 

4.351 

22 

6 

2 

1 

.8202 

.7742 

9.0151 

I 

3.040 

3.402 

^19^1 

5.465 

26 

2 

1 

.7570 

.6633110.9471 

1  1 

3.977 

11.935 
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14 
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.6954 

.6554 

7.576! 

1 

2.557’ 

3.609 

^21^7^1 

1.444 

8 

3 

1 

.7497 

.7232 

5.6821 

1 
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1^/  Following  Strahler's  method 
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Figure  15 

Ellipses  constructed  to  scale  of  four  selected  basins 
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Figure  16 

Comparison  of  the  elliptical  form  of  two  basins  to  the  actual  basin  shapes 
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Comparison  of  the  elliptical  form  of  two  basins  to  the  actual  basin  shapes 


This  value  vjill  be  zero  for  a  circle,  and  approach  one  as  the  ellipse 
approaches  a  straight  line.  Figure  18  shows  the  elliptical  forms  of 
the  four  basins  reduced  to  a  common  minor  axis  so  that  a  direct  com¬ 
parison  may  be  made  of  their  shape. 

Table  13  compares  the  length-width  ratio  and  eccentricity  to  the 
previously  established  elongation  ratio  and  circularity  ratio  for  Spring 
Creek  and  Little  Washita  River  watersheds.  Later,  these  dimensionless 
parameters  will  be  related  to  measured  hydrologic  characteristics  of 
the  basins. 

In  addition  to  the  previously  mentioned  geomorphic  parameters  for 
Spring  Creek,  hypsometric  integrals — the  percent  mass  of  a  basin  above 
the  contour  encompassing  50  percent  of  the  area — have  been  computed. 
Table  14  lists  this  parameter  for  Spring  Creek. 

Table  14. — Hypsometric  integrals  for  basins  of  Spring  Creek  watershed 

_ Basin  identification  number _ Hypsometric  integral _ 

.46 
,62 
.60 
.54 
.54 
,54 


1 

2 

3 

4 

5 

6 


Dr,  T.  A,  Bond  studied  pollen  from  oxbow  deposits  of  the  Washita 
River  as  a  portion  of  his  dissertation  study.  University  of  Oklahoma. 

The  objectives  of  the  investigation  were  to  determine  the  nature, 
successional  history,  environmental  significance  and  relative  ages  of 
oxbow  lake  deposits  in  central  and  southeastern  Oklahoma. 

It  was  concluded  that  spores  and  pollen  are  useful  when  used  in 
conjunction  with  carbon  14  in  dating  sediments  of  oxbow  lakes.  When 
the  pollen  and  spore  assemblages  from  the  oxbow  lakes  were  compared 
with  those  of  the  Domebo  site  it  appeared  that  the  oxbows  are  older 
than  7,000  but  younger  than  10,123  years  B.P,  This  would  place  them 
in  the  late  Wisconsin  age.  An  oxbow  lake  sampled  near  Chickasha  was 
thought  to  be  older  than  either  of  two  oxbow  lakes  sampled  near 
Anadarko,  because  an  assemblage  of  spores  and  pollen  indicative  of  an 
earlier  period  in  the  climatic  succession  were  found  in  the  Chickasha 
oxbov; , 

Observation  wells  have  been  drilled  along  the  valley  in  the  study 
reach  on  1-mile  centers  or  less,  except  for  the  Sugar  Creek  area. 

These  wells  have  provided  enough  information  so  that  an  isopachous  map 
of  river  fill  has  been  completed  except  for  T,  71J,  R.  lOW  (Anadarko), 
and  the  Sugar  Creek  area.  The  volume  of  unconsolidated  sediment  without 
the  two  areas  mentioned  above  is  1.26  cubic  miles  over  an  area  of  127 
square  miles.  T^e  average  thickness  of  fill  over  this  area  is  52.1  feet. 
Once  the  volume  of  sedim-entary  fill  is  knoxTO,  the  volume  of  ground  water 
in  storage  can  be  calculated  monthly  and  compared  directly  with  precipi¬ 
tation  and  other  upstream  variables. 
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Table  13. — Di.tnensionlass 
Little  Washita 

geomorphic  paiame 
watersheds 

ters  for  Spring 
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Basin  ... 
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Circularity  i  Elongation 

.  Length-width 

number 

ratio 

.  ratio 

.  ratio 

Eccentricity 

Spring  Creek 
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.613 

.774 

1.6723 

.8015 

2 

0  0*7 
•  UO  / 

.041 

1.4100 

.  7050 

3 

.777 

.681 

2.6275 

.9247 
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.772 

.636 

2.4745 

.9147 
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.9325 
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Table  13. — Dimensionlsss  geomorphic  parameters  for  Spring  Creek  and 


Littl 

a  Washita  watersheds  (cont’d) 

Basin  . 

identification  j 
number  . 

Circularity 

ratio 

Elongation  ;  L 
ratio  ; 

ength-width  ; 
ratio  ; 

Eccentricity 

Little  Washita 

^11^1 

.7126 

.6193 

2.6052 

.9234 

^12^1 

.3206 

.7919 

1.5947 

.7790 

^13^5^1 

.7916 

.6562 

2.3226 

.9026 

3i4«55i 
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,11  kl 

1.6683 

.0004 
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Figure  18 

Ellipses  of  the  four  selected  basins  reduced  to  a  common  minor  axis  to  compare 
their  shapes 
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CALIBRATION  OF  TRIANGULAR  WEIRS 


Bruce  J.  Blanchard 

Some  measuring  stations  with  vreir  controls  being  used  at  Chickasha 
on  the  subdivided  watershed  studies  are  in  remote  areas  and  are  diffi¬ 
cult  to  reach  in  time  to  calibrate  with  discharge  measurements  in  the 
field.  A  reliable  prediction  equation  for  discharge,  along  with  a 
fev7  discharge  measurements  to  use  as  checks,  would  eliminate  the 
necessity  of  calibration  by  modeling.  The  standard  form  of  the  equation 
for  discharge  over  a  sharp-crested,  V-notch  weir,  Q  =  is  inade¬ 

quate  when  effects  from  variable  approach  velocities,  submergence,  and 
fast-rising  flows  are  encountered,  A  comprehensive  analysis  of  the 
stage-discharge  relationships  of  a  triangular  weir  was  undertaken  to 
improve  the  prediction  of  runoff  from  stream-gaging  stations  with  weir 
controls . 

A  theoretical  approach  to  flow  over  a  sharp-crested,  triangular 
weir  leads  to  the  following  equations 


q  =  ^  Ce  Z 


Vo^  ,  . 
h 


Vo^ 


2g 


Vf 


28 


(1) 


where  q  =  discharge,  =  coefficient  of  contraction,  Z  =  side  slope  of 
the  weir,  a  =  correction  for  nonuniformity  of  the  approach  velocity,  Vq 
mean  approach  velocity,  and  h  =  the  head  causing  flow  over  the  weir. 

This  basic  equation  does  not  take  into  account  the  modification  of 
flow  due  to  submergence,  rate  of  change  of  stage  or  malformation  in 
construction  of  the  weir.  The  effects  of  submergence  on  discharge  over 
a  V-notch  weir  have  been  investigated  in  several  laboratories,  but  all 
experiments  were  on  weirs  less  than  3  feet  wide.  The  weirs  in  use  at 
Chickasha  range  up  to  100  feet  wide,  and  it  was  not  known  if  the  sub¬ 
mergence  relationships  obtained  in  the  laboratory  experiments  vjould 
hold  true. 

Corrections  for  rate  of  change  of  stage  had  been  used  for  some  t-imo 
in  computing  flow  of  large  rivers j  however,  suitable  coefficients  for 
use  on  streams  in  this  area  had  not  been  determined.  Corrections  for 
these  effects  are  made  by  multiplying  the  basic  equation  (1)  by  the 
submergence  factor  equation  (2) ,  and  the  rate  of  change  of  stage 
equation  (3). 

The  effect  of  submergence  on  the  discharge  rate  of  a  weir  was 
assumed  to  be  one  of  superposition  in  x^7hich  actual  flow  xi7as  a  function 
of  the  ratio  of  free  flow  downstream  to  free  flow  upstream  under  the 
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tailwater  head.  The  submergence  factor  T-rhich  is  multiplied  by  the 
free  flow  discharge  has  the  following  form; 


k  1.0 


Cc2  ^2^ 


—  — 
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_  — . 
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Vi2 
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5  ^ 
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“2g 

2^1 

''2g 

(2) 


where  k  is  a  coefficient  usually  found  to  be  l.O,  subscripts  1  and  2 
refer  to  upstream  and  doTnistream  variables  respectively,  and  m  is  a 
constant.  VJhen  h2  is  zero,  the  submergence  factor  is  1.0,  Analysis 
of  the  data  and  the  equations  shows  that  k  should  be  1.0  and  that  a 
could  not  be  predicted  without  plotting  discharge  measurements.  The 
k  was  then  dropped  from  the  equation  and  a  value  of  1.03  x^as  assigned 
to  a  because  it  gave  the  best  fit  of  predictions. 

Rate  of  change  of  stage  over  the  weir  Xiras  assumed  to  modify  the 
flov7  rate  by  the  following  factor 

•^l.O  +  k||  (3) 

dll 

x^hich  is  also  multiplied  times  the  floxi/  and  in  x/hich  is  the  rate  of 
change  of  stage  in  feet  per  minute  and  K  is  considered  constant^ 

In  this  study  the  discharge  data  from  x-7hich  the  final  equations 
were  developed  came  from  an  actual  weir  which  x-jas  slightly  malformed. 

A  correction,  M,  for  this  malformation  xjas  calculated  as  the  ratio  of 
theoretical  discharge  of  a  true-shaped  x-jeir  to  the  theoretical  dis¬ 
charge  of  the  malformed  weir.  The  m.easured  discharge  was  multiplied  by 
this  factor  to  obtain  discharges  for  a  true  x-reir  with  which  a  pre¬ 
diction  equation  could  be  developed.  Therefore,  to  nredict  floX'T  of  any 
actual  xjeir,  the  prediction  equation  must  be  multiplied  by  the 
reciprocal  of  corrections  obtained  in  this  same  manner. 

Coefficients  in  the  above  equations  x-rere  evaluated  by  uaing  82 
discharge  measurements  on  one  vreir.  The  equations  x^^ere  then  fitted 
for  tXAjo  different  conditions;  XArith  and  x^ithout  velocity  data.  The 
resulting  prediction  equations  are; 


^^Corbett,  D.  M. ,  and  others,  1943  Stream  Gaging  Procedure  -  U.  S. 
Geological  Survey  Water  Supply  Paper  888,  page  162,  (Reprinted  1945) 
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Without  velocity  data 


qo  ~  '2.,17>  Z  h]^ 


2.73 


1.0  - 


hi 


2.73 


.405 


1  +  2.13^^^ 


(4) 


dt 


With  discharge  measurement  data 


q  =  2.27  Z  h, 
^o  1 


214 


h  + 


V,2 

i.oa  1 


2g 


2.5 


1.0  - 


h2 


2.71 


hi 


.214 


hi  +  1,08: 


Vi  2 


2g 


2.5 


/ 


1.0  +  .692 


dh(l) 


dt 


M 


(5) 


where  q  =  observed  discharge  in  cfs,  Z  =  least  squares  fitted  weir 
side  slope  ratio  -  horizontal  to  vertical,  h]^  =  upstream  head  in  feet 
above  the  crest  adjusted  for  malformation,  Vi  =  mean  upstream  velocity 
in  feet  per  second,  h2  =  downstream  head  in  feet  above  the  crest  adjust- 

ed  for  malformation,  =  rate  of  change  of  stage  in  feet  per  minute, 

dt 

and  M  =  malformation  coefficient. 

The  two  most  conventional  equations; 

0  =  ch” 

Q  =  ch’^ 

dt 

were  fitted  to  the  data  and  comparison  of  the  standard  errors  indicated 
that  use  of  the  proposed  forms  of  the  discharge  equations  (4)  and  (5) 
are  justified. 


(6) 

(7) 
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ErCPLOMTOriY  STT^DY  OF  THE  REGIME  OF  WaSHITA  RI\'EF  MAIR  STEM  FLOFR 


Dcnn  G.  DeCoursay 

Objectives  of  this  study  x-rere  to  reproduce  a  10-year  regimen  of 
fli^wG  on  the  main  stem  of  the  Washita  River,  then  to  calculata  the 
affects  upon  this  flow  regimen  of  structural  measures  (primarily 
floodvater-retarding  reservoirs)  being  installed  in  the  Washita  River 
Basin  by  the  U.  S.  Department  of  Agriculture,  as  authorized  bv  The 
Flood  Cortrol  Act  of  1944  and  subsequent  amendments'^-.  To  accomplish 
these  objectives,  synthetic  floods  of  surface  runoff  on  tributaries  to 
the  Washita  River  were  developed  then  routed  down  the  main  channel  of 
the  Washita  to  form  the  regime  of  main  stem  flows.  The  unit  hydrograph 
approach,  although  crude,  TTas  selected  as  the  main  tool  for  developing 
the  tributary  inflow.  Records  from  gaging  stations  on  tributaries  to 
the  main  stem  were  used  to  derive  synthetic  unit  hydrographs  for  all 
tributaries  to  the  Washita  River.  Six  of  the  gaged  watersheds  (Barnitz, 
Rush,  Hock,  Little  Washita,  Pond,  and  Sugar)  provided  a  representative 
sample  of  t^e  entrJre  Washita  watershed  and  had  records  suitable  for 
hydrograph  analysis.  Changes  in  the  flew  regime  at  six  main  stem  river 
stations  (Cheyenne,  Clinton,  Carnegie,  Tablet,  Pauls  Valley,  and  Dum^Tood) 
were  observed.  The  period,  1941-50  was  selected  because  data,  was  avail¬ 
able  for  the  entire  period  at  all  main  stem  river  stations  and  there 
were  no  f loodwater-retarding  structures  present.  ‘P'lgure  19  shoTiTs  the 
location  of  the  tributary  and  main  stem  gaging  stations. 

The  following  steps  were  completed  in  reaching  the  ebjeettves; 

1.  Precipitation,  n-noff,  and  topographic  data  were  collected. 


2.  A  mf'dified  form  of  the  Gamma  function  was  used  to  define  the 
■shape  of  the  tributary  unit  hydrographs^ 


3.  The  unit  hydrograph  peak  rate,  time  to  peak,  and  recession 
constant  used  in  the  hydrograph  equation  have  been  related  to  watershed 
characteristics  as  follows? 


tp  =  0.085  W^*^^ 

Pp  =  46.1  tp~“”^^  t,j2.55 

d  =  1.76  In  i)A  -3.85 


1  1-08 
-t-cn 


where  tp  is  the  time  to  peak  in  hours ,  Pp  is  the  peak  rate  in  cf s ,  d  is 
the  recession  constant,  W  is  the  average  watershed  width  in  miles,  l(,g 
is  the  length  of  the  watershed  to  the  point  onposite  the  centroid  in 
miles,  and  DA  is  the  drainage  area  in  square  miles.  The  value  of  tp 
was  rounded  off  to  the  nearest  hour  before  being  used  in  the  prediction 
equation  for  pp. 

this  paper  these  structural  measures  will  hereinafter  be  referred 
to  as  the  FR8  Program. 

1^.\RS  41-116. 
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Figure  19 

Washita  River  sub-watersheds 


) 


4.  The  equation  used  to  calculate  the  tributary  runoff  volume  in 
inches  was: 


Q  =  M  •  P  •  I 


API  +  .068  P 


1.24 


0  is  the  runoff  in  inches,  M  is  a  geologic  constant,  P  is  the  precipi¬ 
tation,  I  is  the  maximum  1-hour  rainfall  intensity,  and  API  is  an 
antecedent  basin  moisture  index  which  is  a  variable  depending  upon 
season.  The  deplctioii  rate  in  the  index  is  a  function  of  the  inverse 
of  the  long-term  temperature  curve  and  varies  continuously  from  0,67 
in  J:’ly-August  to  0.90  in  January, 

Channel  losses  were  assumed  to  var}'  from  reach  to  reach,  but  xjere 
assumed  to  be  constant  within  a  reach.  They  were  taken  as  a  constant 
percent  of  the  flox,7  rate. 

5.  The  Washita  Basin  was  subdivided  Into  70  subwatersheds ,  the 
location  and  physical  characteristics  of  which  are  shown  in  Figure  19 
and  Table  15,  respectively. 


6,  The  successive  average-lag  m.ethod  of  flood  routing  was  used 
to  route  all  unit  hydrographs  including  channel  losses  to  the  main  stem 
stations.  Flood  hydrographs  were  developed  by  sum.ming  the  products  of 
the  routed  unit  hydroyraphs  and  the  precipitation  excess. 


7,  The  effect  of  stri’.ctures  on  local  or  tributary  flov;  was  de¬ 
termined  by  breaking  the  f]ow  Into  two  parts,  uncontrolled  flow  from 
the  drainage  area  below  structures,  and  the  outflow  from  structures. 

The  unit  hydrograph  from  the  uncontrcl3.ed  area  is  a  modification  of  the 
original  unit  hydrograph  for  the  tributary  as  follows: 


t 


u 


1.04 


0,159 


q  =  1,12 
-u 


DAt 


1.33 


where  t,,  and  are  the  times  to  peak  of  the  uncontrolled  area  and  the 
original  tributary  unit  bydrographs  respectively,  q,^^  and  q^.  are  the  pea 

rates  of  the  uncontrolled  area  and  the  original  tributary  unit  hydro¬ 
graphs  respectively,  and  and  DA^  are  the  drainage  areas  of  the 

uncontrolled  area  and  original  tributaries  respectively. 

Using  these  equations,  unit  h^^drographs  for  the  uncontrolled  area 
were  developed  and  routed  to  the  m.ain  stem  stations. 
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Outflox?  from  structures  was  calculated  hy  routing  precipitation, 
runoff,  and  evaporation  through  a  reservoir.  In  order  to  facilitate 
economy  of  operations,  the  physical  characteristics  of  an  average 
structure  of  one-square-mile  drainage  area  were  developed.  These 
characteristics  vxere  based  on  242  structures  in  the  Washita  basin  with 
a  combined  drainage  area  of  800  square  miles. 

Rainfall,  runoff,  and  evaporation  for  the  10-year  period  on  each 
of  the  70  subwatersheds  were  routed  through  the  one-square-mile 
structure  and  the  discharge  multiplied  by  the  number  of  square  miles 
behind  structures.  The  structure  outflow  was  then  routed  down  the 
river  to  the  main  stem  stations. 

Flood  hydrographs  reflecting  the  FRS  program  were  developed  by 
multiplying  the  routed  unit  hydrographs  from  the  uncontrolled  area  by 
pi'cr-lpT  fion  excess  and  adding  the  routed  outflow  from  structures. 

Figure  20  is  a  representative  sample  of  the  flood  hydrographs  with 
and  without  the  FRS  program. 

Discussion 


Flow- duration,  flow-volume,  and  flow-frequency  curves  were 
developed  for  each  river  station  for  the  observed  10-year  period,  the 
predicted  10-year  period  without  structures,  and  the  predicted  10-year 
period  with  structures.  Flow-duration  and  flow-volume  curves  at 
Carnegie  are  presented  in  Figure  21.  The  difference  in  the  observed 
and  predicted  x^rithout  structures  flow-duration  curves  at  low  discharge 
rates  is  due  to  errors  in  separating  base  flow  from  surface  flox<r. 

Since  the  period  of  record  xi;as  short,  a  partial  duration  series 
was  used  to  determine  frequency  curves  from  mean  daily  discharge  rates. 
The  curves  for  Carnegie  are  shown  on  Figure  22.  Because  rate  of  change 
of  the  flow-duration  and  flow-volume  curves  are  more  indicative  of 
changes  that  have  taken  place,  curves  of  this  type  have  been  drawn  for 
each  river  station.  The  curves  for  the  Carnegie  river  station  are 
shown  on  Figures  23  and  24. 

Analysis  of  the  frequency  curves  and  distribution  curves  of  all 
stations  used  in  the  study  have  been  completed.  Follox\ring  are  the 
conclusions  reached: 

Flow-frequency  -  The  period  of  record  used  in  the  study  was  too 
short  to  define  the  changes  in  flood  frequency  caused  by  an  FRS  program 
for  return  periods  greater  than  10  to  15  years . 

Large  storms  are  overpredicted  and  small  storms  underpredicted.  In 
the  upper  reaches  of  the  x^atershed  where  rainfall  is  lighter  and  the 
geologic  constant  lower,  the  predicted  frequency  curve  for  the  10-year 
period  without  structures  plots  lox-xer  than  the  observed  record.  In  the 
lower  reaches,  the  predicted  frequency  curves  for  the  10-year  period 
with  structures  plots  above  the  observed  record.  In  the  central  part 
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Figure  20 

Mean  daily  flow  rates  for  Washita  River  at  Carnegie,  Oklahoma 
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Figure  21 

Flow“duration  and  flow”Vol\’ue  curves  for  Washita  River  at  Carnegie,  Oklahoma 
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of  the  watershed,  the  correspondence  is  good.  The  trend  is  associated 
with  differences  in  rainfall  across  the  watershed  and  with  differences 
in  the  geologic  constants.  These  deviations  will  not  affect  the  rela¬ 
tive  difference  between  the  frequency  curves  of  the  without  and  with 
structures  conditions. 

The  frequency  curve  obtained  from  the  predicted  mean  daily  peaks 
without  structures  is  used  as  a  base  to  calculate  the  percent  reduction 
in  peak  rates  for  each  station  caused  by  the  structures  program.  These 
reductions  for  various  frequencies  of  occurrence  are  shown  in  Table  16. 
Blank  spaces  in  the  table  are  due  to  a  poor  plotting  pattern  of  the 
points.  The  frequency  of  occurrence  was  used  as  the  weight  factor  to 
obtain  the  weighted  average  reduction  for  each  station.  The  percent 
of  the  drainage  basin  that  was  assumed  to  lie  behind  flood  retarding 
iJtructures,  weighted  by  the  volume  of  surface  runoff  at  the  station, 
is  also  shown  in  Table  16.  Figure  25  is  a  plot  of  the  percent  of 
the  watershed  behind  floodwater  retarding  structures  vs.  the  reduction 
in  peak. 

The  1941-50  data  presented  in  Table  16  show  a  relatively  constant 
reduction  in  peaks.  However,  there  are  two  deviations  from  this  trend. 

In  areas  of  low  runoff,  Cheyenne  and  Clinton,  there  is  a  gradual  decrease 
in  the  percent  reduction  in  peaks  with  a  decrease  in  the  frequency  of 
occurrence.  This  is  caused  by  the  com.bination  of  a  low  volume  of  runoff 
and  proportionately  high  evaporation  loss.  The  other  exception  is  in 
the  lox\7er  reaches  at  Pauls  Valley  and  Durwood  where  the  large  events 
show  almost  no  reduction.  This  was  caused  by  emergency  spillway  flow 
which  did  not  occur  at  all  in  the  upstream  reaches  and  only  very  seldom 
in  the  central  portion.  This  would  indicate  a  constant  reduction  in 
peaks  with  respect  to  the  frequency  of  occurrence  unless  the  event  is 
large  enough  to  produce  emergency  spillway  flow. 

Flow-distribution  -  The  distribution  of  flow  duration  by  rate  for 
the  Carnegie  station  is  portrayed  on  Figure  23.  For  this  raid-reach 
station  the  flow-duration  for  the  without  structures  condition  is 
approximately  equally  distributed  about  a  flow  rate  of  70  cfs.  Curves 
for  upstream  stations  (dryer  areas)  show  a  greater  percentage  flow- 
duration  at  smaller  rates,  and  curves  for  dov/nstreara  stations  (wetter 
areas)  show  the  greater  percentage  at  the  larger  rates.  Table  17  shows 
the  mean  annual  surface  flow  rate,  the  discharge  rate  at  which  peak 
duration  occurs,  and  the  ratio  of  the  rate  at  peak  occurrence  to  the 
mean  annual  flow. 

The  distribution  of  flow-duration  by  rate  curves  shows  that  the 
FRS  program  causes  an  increase  in  the  duration  of  flow  in  the  central 
range  of  flow  rates.  There  is  a  corresponding  decrease  in  the  duration 
of  flow  in  the  extreme  upper  range.  At  the  lower  end,  the  durations 
are  approximately  the  same.  The  point  of  intersection,  or  cross-over 
point,  at  about  3,550  cfs,  of  the  two  curves  on  Figure  23  is  a  measure 
of  the  effect  of  structure  installation  on  the  f lox^j-duration  character¬ 
istics  of  the  Washita  River.  Table  17  lists  the  cross-over  point  and 
its  ratio  to  the  mean  annual  flow  for  all  stations.  The  point  at  which 
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Table  16. --Percent  reduction  in  mean  daily  flow  peaks 


Weighted 


Station 

2 

1 

.5 

.2 

.1 

average 

behind  stru( 

Cheyenne 

5C.7 

50.6 

48.3 

43.9 

40.4 

53.9 

62.5 

Clinton 

35.1 

38.9 

47.5 

42.1 

33.7 

33.1 

54.2 

Carnegie 

26.3 

24.6 

24.0 

24.3 

24.3 

25.4 

49.3 

Tabler 

18.9 

18.9 

17.6 

20.0 

21.3 

13.3 

38.1 

Pauls  Valley 

15.4 

24.6 

28.6 

— 

— 

20.8 

41.4 

Durwood 

25.2 

26,6 

24.6 

12.7 

0.0 

24.2 

40.6 

a.J.v:::cL 


*Weighted  by  volume  of  runoff 


Table  17. — Pi stribution  o flow-duration  by  ra t_3 


Station 

Mean 
annual 
surface 
flov7  rata 
(cfs) 

Peak  duration 

Cross-over  point 

Discharge 

rate 

(cfs) 

Ratio  to 

mean 

annual 

Discharge 

rate 

(cfs) 

Ratio  to 

mean 

annual 

Cheyenne 

31 

10 

.32 

1,350 

43.5 

Clinton 

91 

30 

.33 

2,200 

24.2 

Carnegie 

230 

70 

.30 

3,550 

15.4 

Tabler 

435 

70 

.16 

4,500 

10.3 

Pauls  Valley 

645 

300 

.47 

6,300 

9.8 

Durwood 

1,325 

700 

.43 

16,300 

10.0 

•4 
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the  two  curves  separate  at  the  lower  rates  of  flox\T  is  shown  more  clearly 
on  Figure  24.  However,  the  point  at  which  they  separate  is  poorly  de¬ 
fined.  The  ratio  of  the  cross-over  point  to  the  mean  annual  flow  is 
plotted  vs.  the  average  area  (sq.mi.)  inches  of  runoff  in  Figure  26. 

The  FRS  program  causes  an  increase  in  the  volume  of  flow  that 
occurs  at  flow  rates  in  the  central  range,  and  a  decrease  at  the  upper 
extreme  of  flow  rate.  The  distribution  of  flow-volume  by  flow  rate 
for  the  Carnegie  station  is  typical  of  the  Washita  River  and  is 
selected  for  portrayal  on  Figure  24.  Table  18  lists  the  discharge  rate 
at  which  the  peak  volume  of  flow  occurs  for  the  without-structures  con¬ 
dition  and  the  with-structures  condition,  and  the  ratio  of  these  two 
rates.  The  average  of  these  ratios  shows  that  the  FRS  program  X7ill 
reduce  the  rate  at  which  the  peak  volume  occurs  by  about  one- third. 

The  percent  of  flow  that  will  occur  during  this  peak  interval  is 
shown  in  Table  18  for  the  two  watershed  conditions.  The  ratio  of  these 
two  peaks,  also  shown  in  Table  18,  averages  about  1.10.  In  other  words, 
the  volume  of  runoff  during  the  peak  flow  with  the  FRS  program  in  effect 
is  about  10  percent  higher  than  without  the  program  and  occurs  at  a 
lovrer  rate  of  flow. 

The  distribution  of  flow-duration  and  volume  by  rate  curves  shows 
that  intermediate  range  flows  were  more  frequent  and  of  longer  duration 
after  the  FRS  program  was  in  effect.  Changes  resulting  from  these 
structures  at  the  various  river  stations  can  be  obtained  from  Tables 
16-19  and  Figures  25  and  26. 

Volume  of  surface  runoff  -  The  model  used  in  the  study  did  not 
take  into  account  any  seepage  losses  from  the  structures  for  two  reasons 

1.  The  volume-rate  of  seepage  from  such  structures  is  not  known. 

2.  The  destination  of  such  seepage  is  not  known. 

All  other  hydrologic  phenomena  were  taken  into  account. 

Shown  in  Table  19  are  the  average  watershed  inches  of  runoff  per 
year  for  the  10-year  period  without  structures.  Also  shown  In  Table 
19  is  the  ratio  of  the  volume  of  runoff  from  the  watersheds  with  the 
FRS  program  in  effect  to  the  volume  without  a  program.  The 
upstream  reaches  show  about  a  25  percent  reduction  in  total  runoff 
follov/lng  the  construction  of  flood  control  reservoirs,  whereas  the 
downstream  reaches  show  less  than  5  percent.  Explanation  for  the  high 
losses  in  the  upstream  reaches  lies  in  the  coupling  of  high  evaporative 
losses  and  low  surface  runoff  with  a  relatively  large  sediment  pool. 

The  lower  reaches  of  the  basin  have  a  lower  evaporative  loss,  higher 
surface  runoff,  and  relatively  small  sediment  pool. 
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Table  18. --Distribution  of  flow- volume  by  rata 


Station 

Discharge  rate 
at  peak  interval 

Percent  of  flow  occurring 
during  peak  interval 

Untreated 

condition 

(cfs) 

SCS 

condition 

(cfs) 

Ratio 

Untreated 

condition 

SCS 

condition 

Ratio 

Cheyenne 

-- 

-- 

— 

— 

— 

Clinton 

2,500 

1,600 

.  64 

7.5 

8.4 

1.12 

Carnegie 

3,000 

2,500 

.83 

3.3 

8.9 

1.07 

Tabler 

4,700 

3,300 

.70 

7.8 

8.4 

1.03 

Pauls  Valley 

6,700 

4,100 

.61 

7.2 

7.7 

1.07 

Durwood 

19,000 

12,000 

.63 

7.1 

8.1 

1.14 

.68 

1. 10 

Table  19. --Station  runoff 


Station 

Drainage 

area 

(sq.mi.) 

Average 
runoff 
per  year 
(inches) 

Average 
area- inches 
of  runoff 
(sq.mi. in. ) 

Ratio  of  runoff 
treated  to 

untreated  conditions 

Cheyenne 

794 

.52 

413 

.725 

Clinton 

1,977 

.62 

1,226 

.788 

Carnegie 

3,129 

1.00 

3,129 

.880 

Tabler 

4,700 

1.26 

5,930 

.927 

Pauls  Valley 

5,330 

1.65 

3,794 

.951 

Durwood 

7,202 

2.98 

21,462 

.971 
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PERCENT  OF  WATERSHED  BEHIND  STRUCTURES  WEIGHTED  BY  VOLUME  OF  RUNOFF 


Figure  25 

Percent  of  watershed  protected  vs.  reduction  in  peak 
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Figure  26 

Average  area  Inches  of  runoff  vs.  the  ratio  of  the  upper  cross 
over  point  to  the  mean  annual  surface  flow  rate 
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